Probability Functions for the Modulus and Angle of the 
Normal Complex Variate 

By RAY S. HOYT 

This paper deals mainly with various 'distribution functions' and 'cumulative 
distribution functions' pertaining to the modulus and to the angle of the 'normal' 
complex variate, for the case where the mean value of this variate is zero. Also, 
for auxiliary uses chiefly, the distribution function pertaining to the reciprocal 
of the modulus is included. For all of these various probability functions the 
paper derives convenient general formulas, and for four of the functions it supplies 
comprehensive sets of curves; furthur, it gives a table of computed values of the 
cumulative distribution function for the modulus, serving to verify the values 
computed by a different method in an earlier paper by the same author. 1 

Introduction 

IN THE solution of problems relating to alternating current networks 
and transmission systems by means of the usual complex quantity 
method, any deviation of any quantity from its reference value is naturally 
a complex quantity, in general. If, further, the deviation is of a random 
nature and hence is variable in a random sense, then it constitutes a 'complex 
random variable,' or a 'complex variate,' the word 'variate' here meaning 
the same as 'random variable' (or 'chance variable'— though, on the whole, 
'random variable' seems preferable to 'chance variable' and is more widely 
used). 

Although a complex variate may be regarded formally as a single ana- 
lytical entity, denotable by a single letter (as Z), nevertheless it has two 
analytical constituents, or components: for instance, its real and imaginary 
constituents (X and Y); also, its modulus and amplitude (\Z\ and 0). 
Correspondingly, a complex variate can be represented geometrically by 
a single geometrical entity, namely a plane vector, but this, in turn, has 
two geometrical components, or constituents: for instance, its two rec- 
tangular components (X and F); also, its two polar components, radius 
vector and vectorial angle (R — \ Z\ and 6). 

This paper deals mainly with the modulus and the angle of the complex 
variate, 2 which are often of greater theoretical interest and practical im- 

J " Probability Theory and Telephone Transmission Engineering," Bell System Tech- 
nical Journal, January 1933, which will hereafter be referred to merely as the "1933 
paper". 

1 Throughout the paper, I have used the term 'complex variate' for any 2-dimensional 
variate, because of the nature of the contemplated applications indicated in the first 

318 



PROBABILITY FUNCTIONS FOR COMPLEX VARIATE 319 

portance than the real and imaginary constituents. The modulus variate 
and the angle variate, individually and jointly, are of considerable the- 
oretical interest; while the modulus variate is also of very considerable 
practical importance, and the angle variate may conceivably become of 
some practical importance. 

The paper is concerned chiefly with the 'distribution functions' 3 and the 
'cumulative distribution functions' pertaining to the modulus (Sections 3 
and 5) and to the angle (Sections 6 and 7) of the 'normal' complex variate, 
for the case where the mean value of this variate is zero. The distribution 
function for the reciprocal of the modulus is also included (Section 4). 

The term 'probability function' is used in this paper generically to include 
'distribution function' and 'cumulative distribution function.' 

To avoid all except short digressions, some of the derivation work has 
been placed in appendices, of which there are four. These may be found 
of some intrinsic interest, besides facilitating the understanding of the 
paper. 

1. Distribution Function and Cumulative Distribution Function 
in General: Definitions, Terminology, Notation, Relations, 

and Formulas 

The present section constitutes a generic basis for the rest of the paper. 

Let t denote any complex variate, and let p and a denote any pair of 
real quantities determining t and determined by t. (For instance, p and 
a might be the real and imaginary components of t, or they might be the 
modulus and angle of t.) Geometrically, p and a may be pictured as gen- 
eral curvilinear coordinates in a plane, as indicated by Fig. 1.1. 

Let t' denote the unknown value of a random sample consisting of a 
single r-variate, and p' and a' the corresponding unknown values of the 
constituents of r' . 

Further, let G{p, a) denote the 'areal probability density' at any point 
p,a in the p,cr-plane, so that G(p,a)dA gives the probability that r' falls 
in a differential area dA containing the point t; and so that the integral of 

paragraph of the Introduction, and also because the present paper is a sort of sequel to 
my 1933 paper, where the term 'complex variate' (or rather, 'complex chance-variable') 
was used throughout since there it seemed clearly to be the best term, on account of the 
field of applications contemplated and the specific applications given as illustrations. 
However, for wider usage the term 'bivariate' might be preferred because of its prevalence 
in the field of Mathematical Statistics; and therefore the paper should be read with this 
alternative in view. 

3 The term ' distribution function' is used with the same meaning in this paper as in 
my 1933 paper, although there the term ' probability law' was used much more frequently 
than 'distribution function,' but with the same meaning. 
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G(p,o)dA over the entire p,cr-plane is equal to unity, corresponding to 
certainty. 

For the sake of subsequent needs of a formal nature, it will now be as- 
sumed that G{p,a) = at all points p,o outside of the pi , pj , tn , a* quad- 
rilateral region in the p,(r-plane, Fig. 1.1, bounded by arcs of the four heavy 
curves, for which p has the values pi and p 2 and a the values a\ and o-> , 
with pi and v\ regarded, for convenience, as being less than p 2 and cr 2 respec- 
tively. Further, G(p,a) will be assumed to be continuous inside of this 




P+ dp P < 

Pa 

Pi 

Fig. 1.1 — Diagram of general curvilinear coordinates. 

quadrilateral region, and to be non-infinite on its boundary. Hence, for 
probability purposes, it will suffice to deal with the open inequalities 

P, < P < Pi , (1.1) (7i < a < <r 2 , (1.2) 

which pertain to this quadrilateral region excluding its boundary; and thus 
it will not be necessary to deal with the closed inequalities pi ^ p ^ P2 
and <ti ^ a ^ o 2 , which include the boundary. 4 

4 The matters dealt with generically in this paragraph may be illustrated by the fol- 
lowing two important particular cases, which occur further on, namely: 

POLAR COORDINATES: p =\ T \ = R, a = B = angle of t. Then Pl = R x = 0, 
pi = Ri= oo, ai = di = 0, tr 2 = 2 = 2t, whence (1.1) and (1.2) become < R < » 
and < 6 < 2v, respectively. 

RECTANGULAR COORDINATES: p = Re T - *, a = Im r = y. Then p, = *, = 
— °° , P2 = X% = °o, a\ = >'i = — », <r 2 = y» = °°, whence (1.1) and (1.2) become — °o < 
a; < °o and — «>_< y < =», respectively. 
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A generic quadrilateral region contained within the quadrilateral region 
Pi , P2 , <n , 02 in Fig. 1.1 is the one bounded by arcs of the dashed curves 
p 3 , pi , 0-3 , <n , where P3 < Pi and a 3 < 0-4 . Here, as in the preceding 
paragraph, it will evidently suffice to deal with open inequalities. 

Referring to Fig. 1.1, the probability functions with which this paper 
will chiefly deal are certain particular cases of the probability functions 
P(p, a), P(p I (r 34 ) and @(p34 , a 3i ) occurring on the right sides of the follow- 
ing three equations respectively: 

p(p < p' < p + dp, a < a' < a + da) = P(p,a)dpda, (1.3) 

p{p < p' < p + dp, <r 3 < a' < 0-4) = P{p I a 3i )dp, (1.4) 

p{p 3 < P < P\ , ffs < a' < 0-4) = Q(pzi , C34). (1.5) 

These equations serve to define the above-mentioned probability functions 
occurring on the right sides in terms of the probabilities denoted by the 
left sides, each expression p( ) on the left side denoting the probability 
of the pair of inequalities within the parentheses. 5 Inspection of these 
equations shows that: P(p,a) is the 'distribution function' for p and a 
jointly; P(p | 0^4) is a 'distribution function' for p individually, with the 
understanding that a' is restricted to the range 3 -to-o-4 ; Q{pu ,031) is a 
'cumulative distribution function' for p and a jointly. 

Since the left sides of (1.3), (1.4) and (1.5) are necessarily positive, the 
right sides must be also. Hence, as all of the probability functions occur- 
ring in the right sides are of course desired to be positive, the differentials 
dp and da must be taken as positive, if we are to avoid writing | dp \ and 
I da I in place of dp and da respectively. 

Returning to (1.3), it is seen that, stated in words, P(p,a) is such that 
P{p.a)dpda gives the probability that the unknown values p' and a' of 
the constituents of the unknown value r' of a random sample consisting 
of a single r-variate lie respectively in the differential intervals dp and da 
containing the constituent values p and a respectively. Thus, unless 
dpdo is the differential element of area, P(p,a) is not equal to the 'areal 
probability density,' G(p,a), defined in the fourth paragraph of this section. 
In general, if £ is such that EJpda is the differential element of area, then 
p(p ) a) = EG(p, a). (An illustration is afforded incidentally by Appendix A.) 

P{p,a), defined by (1.3), is the basic 'probability function,' in the sense 
that the others can be expressed in terms of it, by integration. Thus 

6 Thus p in p( ) may be read 'probability that' or 'probability of.' 
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P(p | a u ) and P{a \ pu\ denned respectively by (1.4) and by the correlative 
of (1.4), can be expressed as 'single integrals,' as follows 6 : 



7) 



P(fi | (734) = f* P M d °> C 1 - 6 ) P b I " 34) = f 4 PM dp - (L 

Q(pza , ffiOi denned by (1.5), can be expressed as a 'double integral,' funda- 
mentally; but, for purposes of analysis and of evaluation, this will be replaced 
by its two equivalent 'repeated integrals': 

<2W, *«) = j" [jH ?M d «~\ d P = [* [{" p M d p] da > ^- 8) 

the set of integration limits being the same in both repeated integrals 
because these limits are constants, as indicated by Fig. 1.1. On account 
of (1.6) and (1.7) respectively, (1.8) can evidently be written formally 
as two single integrals: 

Q(psi , <ru) = I"' P (P I ™) d P = \ P ( ff I Pa*) d °> ^- 9 ) 

but implicitly these are repeated integrals unless the single integrations in 
(1.6) and (1.7) can be executed, in which case the integrals in (1.9) will 
actually be single integrals, and these will be quite unlike each other in 
form, being integrals with respect to p and a respectively— though of course 
yielding a common expression in case the indicated integrations can be 
executed. 

The particular cases of (1.4) and (1.5) with which this paper will chiefly 
deal are the following three: 

p(p <p' < p + dp, ffl < a' < <r 2 ) = P{p | cm) dp m P (p) dp, (1.10) 

p(pi <p' <p,ai<a' < a 2 ) = (?(< P,<7i 2 ) = Q(p), (1-11) 

Pip < p' < p 2 , *i < a < tf%) = Q(> P,°n) = Q*(p). (1-12) 

"The single-integral formulation in (1.6) can be written down directly by mere inspec- 
tion of the left side of (1.4). Alternatively, (1.6) can be obtained by representing the left 
side of (1.4) by a repeated integral, as follows: 

/p+dp r rat I r r°* 1 

J P{p, a)da \dp - J Pip, <r)d* Up, 

whence (1.6); the last equality in the above chain equation in this footnote evidently 

/x+dx 
f(x)dx = fix)dx, since each side of this equa- 
tion represents dA . the differential element of area under the graph of fix) from x to 
x + dx. 
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In each of these three equations the very abbreviated notation at the ex- 
treme right will be used wherever the function is being dealt with exten- 
sively, as in the various succeeding sections. Such notation will not seem 
unduly abbreviated nor arbitrary if the following considerations are noted: 
In (1.10), an corresponds to the entire effective range of a, so that P(p \ a^) 
is the 'principal' distribution function for p. Similarly, in (1.11), Q(< p,0u) 
is the 'principal' cumultive distribution function for p. In (1.12), the star 
indicates that Q*(p) is the 'complementary' cumulative distribution func- 
tion, since Q(p) + Q*(p) = Q{pn , <ru) = 1, unity being taken as the measure 
of certainty, of course. 

For occasional use in succeeding sections, the defining equations for 
the probability functions pertaining to four other particular cases will 
be set down here: 

p(p<p' <p + dp, in < 1/ < o) = Pip | < a) dp, (1.13) 

p{p<p' < p + dp, a < a' < «r 8 ) = P(p \ > a) dp, (1.14) 

p( Pl <p' <p,a 1 <a' <a) = Q{< p, < a), (1.15) 

Pip < p' < p 2 , <n < a' < a) = Q{> p, < a). (1.16) 

It may be noted that (1.13) and (1.14) are mutually supplementary, in the 
sense that their sum is (1.10). Similarly, (1.15) and (1.16) are mutually 
supplementary, in the sense that their sum is <2(pi?,< o") = (?(< o".Pi2), 
which is the correlative of (1.11). 

This section will be concluded with the following three simple trans- 
formation relations (1.17), (1.18) and (1.19), which will be needed further 
on. They pertain to the probability functions on the right sides of equa- 
tions (1.3), (1.4) and (1.5) respectively, h and k denote any positive real 
constants, the restriction to positive values serving to simplify matters 
without being too restrictive for the needs of this paper. 

P(hp, ka) = 1 P(p,<T), (1.17) 

P{hp\k<r 3i ) =\p(p\ <*»4), (1.18) 

Q(h P3i ,kau) = Q(P34, C734). (1.19) 

Each of the three formulas (1.17), (1.18), (1.19) can be rather easily 
derived in at least two ways that are very different from each other. One 
way depends on probability inequality relations of the sort 

p(l<l'<t+dt) = p{gt<gt'<gt+d[gt]), (1.20) 

p{h<t'<U) = p(gh<gt'<gh), (1.21) 
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where t stands generically for p and for a, and g is any positive real constant, 
standing generically for h and for k; (1.20) and (1.21) are easily seen to be 
true by imagining every variate in the universe of the /-variates to be 
multiplied by g, thereby obtaining a universe of (g/)-variates. A second 
way of deriving each of the three formulas (1.17), (1.18), (1.19) depends on 
general integral relations of the sort 

f /a) it = - r m digo =-[ ob f (~) d\. (i.22) 

Ja g Jga g J oa \g/ 

A third way, which is distantly related to the second way, depends on the 
use of the Jacobian for changing the variables in any double integral; thus, 



P(X,M) 



dXdp. 
dpaa 


= 


d(X,/0 
d(p,cr) 


-1 + 


dip,a) 

d(\ lM ) 



(1.23) 



the first equality in (1.23) depending on the fact that the two sets of vari- 
ables and of differentials have corresponding values and hence are so re- 
lated that 

pip<p'<P+dp, o<o'<<T+d<j) = p{\<\'<\+d\, n<n'<n+dn), (1.24) 



whence 



P(p,a) | dpdtr | = P(X,m) I d\dp |. 



2. The Normal Complex Variate and Its Chief Probability Functions 
The 'normal' complex variate may be defined in various equivalent ways. 
Here, a given complex variate z = x + iy will be defined as being 'normal' 
if it is possible to choose in the plane of the scatter diagram of 2 a pair of 
rectangular axes, u and v, such that the distribution function P(u,v) 
for the given complex variate with respect to these axes can be written in 
the form 8 



P(u,v) = 



1 



exp ■- 



2irouO v [_ ^o 



I ~ 2Sl] 



= P(u)P(v). 



(2.1) 



We shall call w = u + iv the 'modified' complex variate, as it represents 
the value of the given complex variate z = x + iy when the latter is referred 
to the M,7>-axes; P(«) and P{v) are respectively the individual distribution 
functions for the u and v components of the modified complex variate; and 

7 Defined by equation (1.3) on setting p = u and a = v. 

8 This equation is (12) of my 1933 paper. It can be easily verified that the (double) 
integral of (2.1) taken over the entire u, »-plane is equal to unity. 
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S u and S v are distribution parameters called the 'standard deviations' of 
u and v respectively. If / stands for u and for v generically, then 

^.-vbiHrss?]' (2 - 2) * -£'*»* a« 

From the viewpoint of the scatter diagram, the distribution function 
P(u,v) is, in general, equal to the 'areal probability density' at the point 
u,v in the plane of the scatter diagram, so that the probability of falling 
in a differential element of area dA containing the point u,v is equal to 
P(u,v)dA ; similarly, P(a) and P(v) are equal to the component probability 
densities. In particular, the probability density is 'normal' when P(u,v) 
is given by (2.1). 

Geometrically, equation (2.1) evidently represents a surface, the normal 
'probability surface,' situated above the u, f-plane; and P(u, v) is the ordinate 
from any point u,v in the «,t'-plane to the probability surface. 

The M.zi-axes described above will be recognized as being the 'principal 
central axes,' namely that pair of rectangular axs which have their origin 
at the 'center' of the scatter diagram of z = x + iy and hence at the center 
of the scatter diagram of w = u + iv, so that w = 0, and are so oriented, 
in the scatter diagram that uv = (whereas z t^ and xy 7^ 0, in general). 

In equation (2.1), which has been adopted above as the analytical basis 
for defining the 'normal' complex variate, the distribution parameters are 
Sv, and S v ; and they occur symmetrically there, which is evidently natural 
and is desirable for purposes of definition. Henceforth, however, it will be 
preferable to adopt as the distribution parameters the quantities S and b 
defined by the pair of equations 9 

S 2 = Si + Si , (2.4) bS 2 = Si - Si , (2.5) 



whence 

-•2 r.2 < rci let \2 



si + si i + (s v /s u y 



(2.6) 



From (2.4), S is seen to be a sort of 'resultant standard deviation.' The 
last form of (2.6) shows clearly that the total possible range of b is 

— 1 ^ b ^ 1 , corresponding to «> ^ S v /S u ^ 0. 

The pair of simultaneous equations (2.4) and (2.5) give 
2Sl = (1+6)S 2 , (2.7) 2^ = (l-6)5 ? , (2.8) 

which will be used below in deriving (2.11). 

9 Equations (2.4) and (2.6) are respectively (14) and (13) of my 1933 paper. 
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With the purpose of reducing the number of parameters by 1 and of 
dealing with variables that are dimensionless, we shall henceforth deal 
with the 'reduced' modified variate W = U + iV denned by the equation 

W = w/S = u/S + iv/S = U+iV. (2.9) 

Thus we shall be directly concerned with the scatter diagram of W = 
U -\- iV instead of with that of w = u + iv. 

The distribution function P(U,V) for the rectangular components U 
and V of any complex variate W = U + iV is denned by (1.3) on setting 
p = U and a = V; thus, 

P(U,V)dUdV = p{V<U'<U+dU t V<V'<V+dV). (2.10) 

When the given variate s = * + iy W normal, so that the modified variate 
w = u + iv is normal, as represented by (2.1), then, since S is a mere con- 
stant, the reduced modified variate W = U + iF defined by (2.9) will 
evidently be normal also, though of course with a different distribution 
parameter. Its distribution function P(U,V) is found to have the formula 

w - wh* -«{- if* " i^»] = P(P)P(I °' (2 ' n) 

where P(t0 anf l ^0 ) are ^ e component distribution functions: 

i r U 2 1 

. w-^^T««4-t+»J' (2 - 12) 

^ = v^T) exp [-r^-J- - (2 - 13) 

These three distribution functions each contain only one distribution 
parameter, namely b; moreover, the variables U = u/S and V = v'/S are 
dimensionless. 

The distribution function P(R,8) for the polar components R and of 
any complex variate W m R(co& 6 + i sin 6) is defined by (1.3) on setting 
p = R and a = 6; thus 

P(R,e)dRdd = p(R<R'<R-{-dR, 0<0'<0+</0). (2.14) 

For the case where W is 'normal,' it is shown in Appendix A that 

r r — R 2 1 

m ' 6) = xVr=T 2 H i^ 1 * 2 (1 ~ 6 cos 20) J (2 ' 15) 

= V^ exp[-L(l - b cos 20)], (2.16) 

7T 

10 This formula can be obtained from (2.1) by means of (2.7), (2.8), (2.9) and (1.17) 
after specializing (1.17) by the substitutions p = U, a = v and h = k = 1/5. It is (16) 
of my 1933 paper, but was given there without proof. 
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where 

L = R*/{l-b>). (2.17) 

In P(R,6) it will evidently suffice to deal with values of in the first 
quadrant, because of symmetry of the scatter diagram. 

The fact that P(R,6) depends on b as a parameter when W is 'normal' 
may be indicated explicitly by employing the fuller symbol P(R,0;b) 
when desired; thus the former symbol is here an abbreviation for the latter. 

In P(R,8) = P(R, 0; b) it will suffice to deal with only positive values of 
b, that is, with 0^6^ 1 (whereas the total possible range of b is — 1^6^1). 
For (2.15) shows that changing b to —b has the same effect as changing 20 
to tt±20, or to tt/2±0; that is, P(R,9; -o) = P(R, tt/2±0; o). 

Seven formulas which will find considerable use subsequently are obtain- 
able from the integrals corresponding to equations (1.13) to (1.16), by setting 
p = R and a = or else p = and a = R, whichever is appropriate, and 
thereafter substituting for P(R,6) the expression given by (2.16), and 
lastly executing the indicated integrations wherever they appear possible. 11 
The resulting formulas are as follows: 



P(R | < 0) = X± exp(-L) ( exp(bL cos 20) dd, 

IT JO 



(2.18) 



P(, | < R) - VT^ l-exp[-L ( l-6cos20)] 
27r 1—6 cos 20 



Pi e | > jQ = viz- -PL-mi - . cos an _ 

27r l—o cos 20 



\ > R) = V 1 ~ ° 2 exp[-£(l - b cos J0)] 
' 2r 1 - 6 cos 20 

Q« R, < 0) = i J" VI exp(-L) / exp(oZ, cos 20) dd\ dR (2.21) 



Vl - 6 2 r tf 1 - exp[-£(l - 6 cos 20)] 



ri-expt-y-t^M)) 

Jo 1 — 6 cos 20 ' v ' 



2tt Jo l-o cos 20 



Q(> /?, < 0) = - f \ VI exp(-Z) / exp(6L cos 20) <# |<«2 (2.23) 

VF^l* r" eX p(oL cos 20) Jn ..... 

= — o exp(-L) / -^— —Ldd. (2.24) 

27T Jo 1-0 cos 20 

Formulas (2. 21) to (2.24) are obtainable also by substituting (2.18) to 
(2.20) into the appropriate particular forms of (1.9). 
When a 0-range of integration is 0-to-^(7r/2), where q = 1, 2, 3 or 4, this 

11 Except that in (2.22) the part 1/(1 - b cos 20) is integrable, as found in Sec. 7, 
equations (7.6) and (7.7). 
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range can be reduced to 0-to-ir/2 provided the resulting integral is mul- 
tiplied by q; that is, 

r a(W2) r rli 

F(d)dd = q F(d)dd, (2.25) 

Jo •'o 

because of symmetry of the scatter diagram. 

3. The Distribution Function for the Modulus 

The distribution function P(R | ]2 ) = P(R) for the modulus R of any 
complex variate W = R(cosd + i sin 0) is denned by equation (1.10) on 
setting p = R, a = 0, <n = 0i = and <r 2 = 02 = 2tt; thus 

P(R)dR = p(R<R'<R+dR, O<0'<2tt). (3.1) 

An integral formula for P(R) is immediately obtainable from (1.6) by 
setting p = R, a = 6, a 3 = oi = 8i = and <r 4 = era = 02 = 27r; thus 

P(J?) - / P(R,e) do. (3.2) 

Jo 

The rest of this section deals with the case where IV m R(cos + i sin 0) 
is 'normal.' Since this case depends on b as a parameter, P(R) is here an 
abbreviation for P(R; b). A formula for P(R; b) can be obtained by sub- 
stituting P(R, 0) from (2.15) into (3.2) and executing the indicated integra- 
tion by means of the known Bessel function formula 

/ eitp(77 cos \{/) d\f/ = irIo(rj), (3.3) 

Jo 

/o( ) being the so-called 'modified Bessel function of the first kind,' of 
order zero. 12 The resulting formula is found to be 13 

This can also be obtained as a particular case of the more general formula 
(2.18) by setting = 2t in the upper limit of integration and then apply- 
ing (3.3). 

In P(R;b) it will suffice to deal with positive values of b, that is, with 
Ogigl, as (3.4) shows that P(R; -b) - P(R;b). 

12 It may be recalled that I (z) = /o(«), and in general that I n (z) = i- n J n (iz). 

In the list of references on Bessel functions, on the last page of this paper, the modified 
Bessel function' is tfeated in Ref. 2, p. 20; Ref. 3, p. 102; Ref. 4, p. 41; Ref. 1, p. 77. 

Regarding formula (3.3), see Ref. 1, p. 181, Eq. (4), v = 0; Ref. 1, p. 19, Eq. (9), fourth 
expression, v = 0; Ref. 2, p. 46, Eq. (10), « = 0; Ref. 3, p. 164, Eq. 103, n = 0. 

"This formula was given in its cumulative forms, / P(R; b)dR, as formulas (.Sl-A) 

and (53-A) of the unpublished Appendix A to my 1933 paper. 
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It will often be advantageous to express P R -,b in terms of b and one or 
the other of the auxiliary variables L and T defined by the equations 

L = rh" (3 - 5) r = Ji = CT w 

Formula (3.4) thereby becomes, respectively, 

P(R;b) = 2-v/Z exp(-L)I Q (bL), (3.7) 

P(R;b) m 2 Jt exp[^J h(T). (3.8) 

Formula (3.8) will often be preferable to (3.7) because the argument of 
the Bessel function in (3.8) is a single quantity, T. 

Because tables of h(X) are much less easily interpolated than tables of 
Mq(X) defined by the equation 

Mo(X) = exp(-X)I (X), (3.9) 

extensive tables of which have been published, 14 it is natural, at least for 
computational purposes, to write (3.4) in the form 

2R r -r 2 i r ^ i 

For use in equation (3.16), it is converient to define here a function 
M\{X) by the equation 

M X {X) = exp(-.Y)/i(X), (3.11) 

corresponding to (3.9) defining M {X). Mi(X) has the similar property 
that it is much more easily interpolated than is Ii(X); and extensive tables 
of Mi(X) are constituent parts of the tables in Ref. 1 and Ref. 6. 

The quantity bR-/{\ — b 2 ) = T, which occurs in (3.4) and (3.8) as the 
argument of 7o( ), and in (3.10) as the argument of M { ), evidently 
ranges from to «> when R ranges from to °o and also when b ranges 
from to 1. Formula (3.10) is suitable for computational purposes for all 
values of the above-mentioned argument bR-Kl — b 7 ) = T not exceeding 
the largest values of X in the above-cited tables in Ref. 1 and Ref. 6. For 
larger values of the argument, and partiularly for dealing with the limiting 

" Ref. 1, Table II (p. 698-713), for X = to 16 by .02. Ref. 6, Table VIII (p. 272- 
283), for X = 5 to 10 by .01, and 10 to 20 by 0.1. Each of these references conveniently 
includes a table of exp(.Y) whereby values of 7o(.V) can be readily and accurately evalu- 
ated if desired. Values of 7o(A") so obtained would enable formulas (3.4), (3.7) and (3.8) 
of the present paper to be used with high accuracy without any difficult interpolations, 
since the table of exp(X) is easily interpolated by utilizing the identity exp(Xi + X 2 ) = 
exp(Xi) exp(A' 2 ). 
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case where the argument becomes infinite, formula (310)— and hence (3.4)— 
may be advantageously written in the form 

2 V -R 2 ~\ T bR 2 1 , s 

m > b) = ^m Hjt»J No |r=?J ■ (3 - 12) 

where 

A 7 (X) - -v/2^Xexp(-X)/ (X) = \^|M.(X), (3.13) 

an extensive table of which has been published. 16 The natural suitability 
of the function No(X) for dealing with large values of X is evident from 
the structure of the asymptotic series for N (X), for sufficiently large values 
of X, which runs as follows: 16 

l 2 l 2 3 a l*tf& 

No(x) ~ 1 + mx + 2m)> + mx)>+'''' (3 - 14) 

whence it is evident that 

JVo(«0 = 1. (3.15) 

For use in Appendix C, it is convenient to define here a function Ni(X) 
by the equation 17 



Nl (X) = \/2irXexp(-X)h(X) = y/lvXM^X), (3.16) 

corresponding to (3.13) defining N (X), with Mi(X) defined by (3.11). 
The asymptotic series for Ni{X), which will be needed in Appendix C, is 18 

whence it is evident that 

tfi(co) * 1. (3.18) 

When b is very nearly but not exactly equal to unity, so that 



bR 2 R 2 R 2 



(3.19) 



1 _ b 2 ~ i - b 2 2(1 - b) ' 

it is seen from (3.4) that P{R;b) is, to a very close approximation, a function 

» Ref. 7, pp. 45-72, for X = 10 to 50 by 0.1, 50 to 200 by 1, 200 to 1000 by 10, and 
for various larger values of X. 

18 Ref. 1, p. 203, with {v, m) denned on p. 198; Ref. 5, p. 366; Ref. 2, p. 58; Ref. 3, p. 
163, Eq. 84; Ref. 4, pp. 48, 84. 

17 N t (X) is tabulated along with No(X) in Ref. 7 already cited in connection with equa- 
tion (3.13). 

18 Ref. 1, p. 203, with (v, m) denned on p. 198; Ref. 5, p. 366; Ref. 2, p. 58; Ref. 3, p. 
163, Eq. 84. 
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of only a single quantity, which may be any one of the three very nearly 

equal expressions in (3.19)— but the last of them is evidently the simplest. 

Fig. 3.1 gives curves of P(R;b), with the variable R ranging continuously 




Fig. 3.1 — Distribution function for the modulus (R = to 2.8). 

from to 2.8 and the parameter b ranging by steps from to 1 inclusive, 
which is the complete range of positive b. Fig. 3.2 gives a*n enlargement 
(along the i?-axis) of the portion of Fig. 3.1 between R = and R *= 0.4, 
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DISTRIBUTION FUNCTION, P(R;b) 

Fig. 3.2— Distribution function for the modulus (22 = to 0.4). 
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and includes therein curves for a considerable number of additional values 
of 6 between 0.9 and 1 so chosen as to show clearly how, with 6 increasing 
toward 1, the curves approach the curve for b = 1 as a limiting particular' 
curve; or, conversely, how the curve for b = 1 constitutes a limiting par- 
ticular curve — which, incidentally, will be found to be a natural and con- 
venient reference curve. This curve, for b = 1, will be considered more 
fully a little further on, because it is a limiting particular curve and be- 
cause of its resulting peculiarity at R = 0, the curve for b = 1 having at 
R = a projection, or spur, situated in the P{R;b) axis and extending from 
0.7979 to 0.9376 therein (as shown a little further on). 

The formulas and curves for b = and 6=1, being of especial interest 
and importance, will be considered before the remaining curves of the set. 

For the case 6 = 0, formula (3.4) evidently reduces immediately to 

P(R;0) = 2R exp (-T? 2 ). (3.20) 

This case, 6 = 0, is that degenerate particular case in which the equiprob- 
ability curves in the scatter diagram of the complex variate, instead of 
being ellipses (concentric), are merely circles, as noted in my 1933 paper, 
near the bottom of p. 44 thereof (p. 10 of reprint). 

For the case 6=1, the formula for the entire curve of P(R; 6) = P(R;1), 
except only the part at R = 0, can be obtained by merely setting 6=1 
in 19 (3.12) as this, on account of (3.15), thereby reduces immediately to 



2 



(-?]■ 



P'(R; 1) = ^/g exp| - , I , (R * 0), (3.21 

P'(R;\) denoting the value of P(R;b) when 6 = 1 but R ^ 0, the restriction 
R 7^ being necessary because the quantity i? 2 /(l — 6 2 ) in (3.12) — and in 
(3.4) — does not have a definite value when 6 = 1 if R = 0. Thus, in Figs. 
3.1 and 3.2, the curve of P'(R;1) is that part of the curve for 6 = 1 which 
does not include any point in the P(R; 6) axis (where R = 0) but extends 
rightward from that axis toward R = +00. The curve of P'(R;l) is the 
'effective' part of the curve of P(R;\), in the sense that the area under the 
former is equal to that under the latter, since the part of the curve of 
P(R;l) at R = can have no area under it. 

P(0;1) denoting (by convention) the value, or values, of P{R;b) when 
R = and 6=1, that is, the value, or values, of P(R;1) when R = 0, it 
is seen, from consideration of the curves of P(R;b) in Figs. 3.1 and 3.2 when 
6 approaches 1 and ultimately becomes equal to 1, that the curve of P(0;1) 
consists of all points in the vertical straight line segment extending upward 
in the P(R;b) axis, from the origin to a height 0.9376 [= Max P(i?;l)], 20 

19 Use of (3.12) instead of (3.4), which is transformable into (3.12), avoids the indefinite 
expression °° .0. «> which would result directly from setting b = 1 in (3.4). 

20 As shown near the end of Appendix B, Max P(R;l) is situated at 22 = and is 
equal to 0.9376. 
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together with all points in the straight line segment extending downward 
from the point at 0.9376 to the point at 0.7979 [= 2/ y/lir = P'(R;D for 
r - o+j. The curve of P(0; 1), because it has no area under it, is the 
'non-effective' part of the curve of P(R;l)- 

Starting at the origin of coordinates, where R = 0, the complete curve 
of P(R;1) consists of the curve of P(0;1), described in the preceding para- 
graph, in sequence with the curve of P'(R;i), given by (3.21). Thus the 
' complete curve of P(R;\) is the locus of a tracing point moving as follows: 
Starting at the origin of coordinates, the tracing point first ascends in the 
P(R; b) axis to a height 0.9376 [= Mux P(R;\)]; second, descends from 
0.9376 to 0.7979 [= 2/ \/2ir = P'(R;D for R = 0+]; and, third, moves 
rightward along the graph of P'(R;1) [b = 1] toward R = + cc . The locus 
of all of the points thus traversed by the tracing point is the complete 
curve 21 of P(R;l). 

In addition to being the principal part ('effective' part) of the curve of 
P(R;1), the curve of P'(R;l), whose formula is (3.21), has a further impor- 
tant significance. For the right side of (3.21), except for the factor 2, will 
be recognized as being the expression for the well-known 1-dimensional 
'normal' law; the presence of the factor 2 is accounted for by the fact that 
the variable R = | i? | can have only posiive values and yet the area under 
the curve must be equal to unity. This case, b = 1, is that degenerate 
particular case in which the equiprobability curves, instead of being ellipses, 
are superposed straight line segments, so that the resulting 'probability 
density' is not constant but varies in accordance with the 1-dimensional 
'normal' law (for real variates), as noted in my 1933 paper, at the top of p. 45 
thereof (p. 11 of reprint). 

All of the curves of P(R;b), where 0^b^\, pass through the origin, 
the curve of P(R;l) [b = 1] being no exception, since the part P(0;1) passes 
through the origin. 

Formula (3.12), supplemented by (3.15), shows that P(R; b) = at 
R = oo ; and this is in accord with the consideration that the total area 
under the curve of P(R;b) must be finite (equal to unity). 

Since P{R\b) = at R = and at R = « , every curve of P{R;b) must 
have a maximum value situated somewhere between R = and R = <» — 
as confirmed by Figs. 3.1 and 3.2. These figures show that when b increases 
from to 1 the maximum value increases throughout but the value of R 
where it is located decreases throughout. 

The maxima of the function P(R;b) and of its curves (Figs. 3.1 and 3.2) 
are of considerable theoretical interest and of some practical importance. 

21 The presence, in the curve of P(R; 1), of the vertical projection, or spur, situated in 
the P{R; b) axis and extending from 0.7979 to 0.9376 therein, is somewhat remindful 
(qualitatively) of the ' Gibbs phenomenon' in the representation of discontinuous periodic 
functions by Fourier series. 
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The cases b = and b = 1 will be dealt with first, and then the general 
case (b = b). 

For the case b = it is easily found by differentiating (3.20) that P(R;b) = 
P(R;0) is a maximum at R = 1/ \/2 = 0.7071 and hence that its maximum 
value is \/2 exp (—1/2) = 0.8578, agreeing with the curve for b = in 
Fig. 3.1. 

For the case 6=1, which is a limiting particular case, the maximum 
value of P(R;b) = P(R;l) apparently cannot be found driectly and simply, 
as will be realized from the preceding discussion of this case. Near the 
end of Appendix B, it is shown that the maximum value of P(R;l) occurs at 
R — (as would be expected) and is equal to 0.9376. This is the maximum 
value of the part P(0;1 of P(R;l). The remaining part of P(R;l), namely 
P'(R;l), whose formula is (3.21), is seen from direct inspection of that 
formula to have a right-hand maximum value at R = 0+, whence this 
maximum value is 2/\/27r = 0.7979. 

For the general case when b has any fixed value within its possible positive 
range (0^6^ 1), it is apparently not possible to obtain an explicit expression 
(in closed form) either for the value of R at which P(R;b) has its maximum 
value or for the maximum value of P(R;b); and hence it is not possible to 
make explicit computations of these quantities for use in plotting curves of 
them, versus b, of which they will evidently be functions. However, as 
shown in Appendix B, these desired curves can be exactly computed, in an 
indirect manner, by temporarily taking b as the dependent variable and 
taking T, defined by (3.6), as an intermediate independent variable. For 
let Re denote the critical value of R, that is, the value of R at which P(R;b) 
has its maximum value; and let Tc denote the corresponding value of T, 
whence, by (3.6), 

T c = bR%/(l-b 2 ). (3.22) 
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Fig. 3.3 — Functions relating to the maxima of the distribution function for the modulus. 
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Then, computed by means of the formulas derived in Appendix B, Fig. 3.3 
gives a curve of Rc and a curve of Max P(R;b), each versus b. Since the 
curve of R c cannot be read accurately at b « 1, there is included also a 
curve of Rc/^/l — b 2 , from which R c can be accurately and ea sily com- 
puted for any value of b \ inci dentally, the curve of R c /\/l — b 2 is simul- 
taneously a curve of VT c /b, on account of (3.22). From Fig. 3.3 it is 
seen that R c varies greatly with b but that Max P R - b varies only a little, 
as also is seen from inspection of Figs. 3.1 and 3.2 giving curves of P(R;b) 
as function of R with b as parameter. 

In Fig. 3.3, the curve of R c shows that for b = 1 the maximum of P(R;b) 
occurs at R = 0; and the curve of Max P{R;b) shows that Max P(R;l) « 
0.94, agreeing to two significant figures with the value 0.9376 found near 
the end of Appendix B. 

4. The Distribution Function for the Reciprocal of the Modulus 

At first, let R denote any real variate, and P(R) its distribution function. 
Also let r denote the reciprocal of R, so that r = 1/R; and let P{r) denote 
the distribution function for r. Then - 

P(r) = R 2 P(R) = P{R)/r\ (4.1) 

If P(R) depends on any parameters, P(r) will evidently depend on the 
same parameters. 

The rest of this section deals with the case where W = R(cos + i sin 0) 
is 'normal.' Since this case depends on b as a parameter, P(R) and P(r) 
are here abbreviations for P{R;b) and P(r;b) respectively. 

As P(R;b) has the distribution function given by (3.4), the distribution 
function for r will be 

p(r ' b) - wr=w exp [o^v] /o [ (T=lv] ' (4 - 2) 

obtained from the right side of (3.4) by changing R to 1/r and multiplying 

22 For if r and R denote any two real variates that are functionally related, say F{r, R) 
= 0, and if dr and dR are corresponding small increments, then evidently 



P(r) 
P(r) \dr\ - P(R) | dR | whence — 



dR 

dr 



dF/dr 



dF/dR 



In particular, if r = 1/R, whence F — r— 1/R, then (4.1) results immediately. 

For a somewhat different and more detailed treatment of change of the variable in 
distribution functions, see Thorton C. Fry, "Probability and its Engineering Uses," 
1928, pp. 153-155. (Cases of more than one variate are treated on pp. 155-174 of the 
same reference.) 
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the result by 1/r 2 , in accordance with (4.1). Evidently P(r;—b) 
= P{r;b). 
By means of (4.1), formulas (3.7) and (3.8) give, respectively, 

P(r;b) = 2(l-V)L 3 ' 2 exp(-L)I (bL), (4.3) 

P(r;b) = 2(1 - b 2 ) [J]'* exp[^] I (T), (4.4) 

wherein L and T are defined by (3.5) and (3.6) respectively, but will now 
be written in the equivalent forms 

which are evidently more suitable for the present section. 

A few particular cases that are especially important will be dealt with 
in the following brief paragraph, ending with equation (4.8). 

For the two extreme values of r, namely and =c , P(r;b) is zero for all 
values of b in the b- range (O^i^l). 

When b = 0, 

P(r;b) = P(r;0) « 1 aq^Z. 1 ]. (4.7) 

P M )=P (r;1) = -^^exp[^]. (4.8) 

Fig. 4.1 gives curves of P(r;b), with the variable r ranging continuously 
from to 1.4 and the parameter b ranging by steps from to 1; however, 
in the r-range where r is less than about 0.6, alternate curves had to be 
omitted to avoid undue crowding. Fig. 4.2 gives an enlargement of the 
section betwen r = 0.2 and r = 0.5, and includes therein the curves that 
had to be omitted from Fig. 4.1. 

In Fig. 4.1 it will be noted that with the scale there used for P(r;b) the 
values of P(r;b) are too small to be even detectable for values of r less 
than about 0.25. Even in the enlargement supplied by Fig. 4.2, the values 
of P(r;b) are not detectable for r less than about 0.2. 

The curves of P{r\b) in Figs. 4.1 and 4.2 would have had to be computed 
from the lengthy formula (4.2) — or its equivalents — except for the fact 
that curves of P(R;b) had already been computed in the preceding section 
of the paper. The last circumstance enabled the P(r;b) curves to be 
obtained from the P(R;b) curves by means of the very simple relation (4.1). 

It will be observed that each curve of P(r;b) [Fig. 4.1] has a maximum 



When 6 = 1, 
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ordinate, whose value and location depend on b. When b increases from 
to 1, the maximum ordinate decreases throughout but the value of r where 
it is located remains nearly constant, at about 0.82, until b becomes about 



0,85 




0.2 0.3 



0.4 0.5 0.6 0.7 0.8 0.9 1.0 

RECIPROCAL OF THE MODULUS, T 



Fig. 4.1— Distribution function for the reciprocal of the modulus (r = to 1.4). 



0.7, after which the location of the maximum value moves rather rapidly 
to about 0.71 fori = 1. 

For the cases b = and b = 1, it is easily found, by differentiating (4.7) 
and (4.8), that the maximum ordinates are located at r = \/2/3 = 0.8165 
and at r = l/\/2 = 0.7071 respectively; and hence, by (4.7) and (4.8), 
that the values of these maximum ordinates are (3-^/3/2 exp(— 3/2) = 
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0.8198 and (4/V2ir) exp (-1) = 0.5871 respectively. These results for 
the cases 6 = and b = 1 agree with the corresponding curves in Fig. 4.1. 



0.40 




0.20 



0.24 



0.28 0.32 0.36 0.40 

RECIPROCAL CF THE MODULUS, r 



0.44 



0.48 



Fig. 4.2 — Distribution function for the reciprocal of the modulus (r = 0.2 to 0.5). 



For the general case where b has any fixed value in the 6-range (0^6^ 1), 
it is apparently not possible to obtain an explicit expression (in closed form) 
either for the value of r at which P(r;b) has its maximum value or for the 
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maximum value of P(r\b). However, as shown in Appendix C, curves of 
these quantities versus b can be computed, in an indirect manner, by 
temporarily taking b as the dependent variable and taking T, denned by 
(4.6), as an intermediate independent variable. For let r c denote the 
critical value of r, that is, the value of r at which P(r;b) has its maximum 
value; and let T c denote the corresponding value of T, whence, by (4.6), 



T c = b/{\-V)r\. 



(4.9) 



Then, computed by means of the formulas derived in Appendix C, Fig. 4.3 
gives a curve of r c and a curve of Max P{r;b), each versus b. From these 
curves it is seen that r c and Max P(r;b) do not vary greatly with b, as also 
is seen from inspection of Fig. 4.1 giving curves of P(r;b) as function of r 
with b as parameter. 
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Fig. 4.3— Functions relating to the maxima of the distribution function for the reciprocal 

of the modulus. 

5. The Cumulative Distribution Function for the Modulus 

The cumulative distribution function Q(<R,6n) ■ Q(R) for the 
modulus R of any complex variate W m R{cos 8 + i sin 0) is denned by 
equation (1.11) on setting p = R, a = 6, Pl = R\ = 0, tn = 0i = and 
<T 2 = 0s = 2t; thus 

Q(R) = p(0<R'<R,0<d'<2T). (5.1) 

Similarly, from (1.12), the complementary cumulative distribution function 
Q(>R,du) = Q*(R) is denned by the equation 



Q*(R) - p(R<R'<™,0<6'<2ir). 



(5.2) 



Q*(R) is usually more convenient than Q(R) for use in engineering ap- 
plications, because it is usually more convenient to deal with the relatively 
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small probability of exceeding a preassigned rather large value of R than to 
deal with the corresponding rather large probability (nearly equal to 
unity) of being less than the preassigned value of R. 

A 'double integral' for Q(R), in the form of two 'repeated integrals,' 
can be written down directly by inspection of the p{ ) expression in 
(5.1) or by specialization of (1.8); thus 

Q(R) = l\l p ( R >°) de \ dR = I [jf p ( R & dR ~\ de - ( s - 3 ) 

Evidently these can be written formally as two 'single integrals,' 

Q(R) - / P(R) dR = [ P{6 I < R) dd, (5.4) 

Jo Jo 

by means of the distribution functions P(R) = P(R | 12 ) and P(6\ <R) 
given by the formulas 

P(R) = [ P(R,d)dd, (5.5) P(fi\ < R) = [ P(R,6)dR. (5.6) 
Jo Jo 

(5.5) is the same as (3.2). (5.6) is a special case of (i.6), and the left side 
of (5.6) is a special case of P(p \ <a) denned by (1.13). 

Similarly, from (5.2), we arrive at the following formulas corresponding 
to (5.3), (5.4), (5.5), and (5.6) respectively: 

Q*(R) = f^ U P(R,d) dd 1 dR = I" [7 P(R,0) dR~\d0, (5.7) 

Q*(R) = [ P(R) dR= I P(6\> R) dd, (5.8) 

Jh Jo 

P(R) = [ P(R,6)dd, (5.9) P(6\ > R) = f P(R,d) dR. (5.10) 

JO J K 

The rest of this section deals with the case where W = R(cos 6 + i sin 6) 
is 'normal.' 23 Since this case depends on b as a parameter, Q{R) and Q*(R) 
are here abbreviations for Q(R;b) and Q*(R;b) respectively. 

A natural and convenient way for deriving formulas for Q(R) is afforded 
by the general formula (5.4) together with the auxiliary general formulas 
(5.5) and (5.6), beginning with the two latter. 

For the 'normal' case, P(R,9) is given by (2.15). When this is sub- 
stituted into (5.5) and (5.6), it is found that each of the indicated integra- 

23 For the 'normal' case, the cumulative distribution function was treated in a very 
different manner in my 1933 paper and its unpublished Appendix A. That paper included 
applications to two important practical problems, and its unpublished Appendix C treated 
a third such problem. (The unpublished appendices, A, B and C, are mentioned in foot- 
note 3 of the 1933 paper.) 
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tions can be executed, giving the two previously obtained formulas (3.4) 
and (2.19) for P(R) m P(R;b) and P(d\ <R) respectively. When these 
are substituted into (5.4), there result two types of single-integral formulas 
for Q(R): A primary type, involving an indicated integration as to R; and 
a secondary type, involving an indicated integration as to 0. Formulas 
of these two types for Q(R) will now be derived. 

An integral formula of the primary type for Q(R) = Q{R;b) can be ob- 
tained by substituting P(R) = P(R;b) from (3.4) into the first integral in 
(5.4), giving 

ot = 2 f vr^ exp [r=^ 2 ] h [r=y *■ (5-11) 

This can also be obtained as a particular case of the more general formula 
(2.21) by setting 6 = 2r in the upper limit of integration and then apply- 
ing (3.3). m .„.,!* 

In (5.11), X is used instead of R as the integration variable in order to 
avoid any possible confusion with R as an integration limit. Thus the 
integrand is a function of X with b as a parameter. Evidently Q{R;b) = 
Q(R;-b). Formula (5.11) is evidently suitable for evaluation of Q(R) by 
numerical integration. 24 

By suitably changing the variable in (5.11), we arrive at the following 
various additional formulas, which, though equivalent to (5.11), are very 
different as regards the integrand and the limits of integration. As previ- 
ously, L denotes R 2 /(\-b 2 ). 

w = ^ f 4f^J 7 ° [i^p] ""■ (5 - 12) 

Q(R) = VT^~b 2 f exp(-X) h(b\)d\, (5.13) 

Q(R) = L\/T^T 2 [ exp(-LX) h(bL\) d\, (5.14) 

Jo 

Q(R) = VT^~b°- [ Ub log X) d\. (5.15) 

These four additional formulas are of some theoretical interest, but ap- 
parently they are less suitable than (5.11) for numerical integration with 
respect to R. A formul a differ ing slightly from (5.11) could evide ntly be 
obtained by taking X/Vl - b 2 as a new variable, and hence R/Vl - b 2 
as the upper limit of integration. 

Corresponding formulas for Q*(R) = Q*(R;b) can of course be obtained 
from the preceding formulas (5.11) to (5.15) inclusive for Q(R) = Q{R;b) 

u In this connection, Appendix D may be of interest. 
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« 
by merely changing the integration limits correspondingly — for instance, 
in (5.11), from 0, R to R, » ; in (5.13), from 0, L to L, « ; and so on. How- 
ever, the first four formulas for Q*(R) so obtained would suffer the disad- 
vantage of each having an infinite limit of integration, rendering those 
formulas unsatisfactory for numerical integration purposes. This difficulty 
can be avoided by making the substitution R = l/r in each of those formulas 
for Q*{R). The resulting formulas are the following five, corresponding to 
(5.11) to (5.15) respectively: 24 

<*•« = ^^ f h exp [ - 3 u GO *- (5 - i8) 

Q*(R) = Ly/T=~V j[ ' 1 exp[ - |J J„ [y] A, (5.19) 

i»exp(— t) 

<2*(7c) = Vl - 6 2 / /o(6 log X) tfX. (5.20) 

Jo 



As a check on (5.16), it is obtainable from (4.2) by integrating the latter 
as to r. 

For purposes of evaluation by numerical integration, formulas (5.11) 
to (5.15) inclusive may evidently differ greatly as regards the amount of 
labor involved and the numerical precision practically attainable. In 
each of these formulas except (5.14) the integrand contains only one param- 
eter, b, while the integration range involves either R or L = R-/(\ — b-). 
In (5.14) the integrand contains two independent parameters, b and L, 
while the integration range is a mere constant, 0-to-l. Similar statements 
apply to formulas (5.16) to (5.20) inclusive. 

A partial check on any formula for Q(R) can be applied by setting R = co , 
since Q(°c) should be equal to unity (representing certainty). If, for 
instance, this procedure is applied to formula (5.13), the right side is lound 
to reduce to unity by aid of the known relation 

[ exp (-A\) J (B\) d\ = . 1 (5.21) 

•'o \/A 2 4- B 2 

together with h(B^) = J (iB\). 

An integral formula of the secondary type for Q*(R) = Q*(R;b) can be 
obtained by substituting (2.20) into the last integral in (5.8), utilizing (2.25), 

» Ref. 1, p. 384, Eq. (1); Ref. 2, p. 65, Eq. (2); Ref. 4, p. 58, Eq. (4.5). 
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changing the variable of integration by the substitution 6 = 0/2, and 
rearranging; thus it is found that 

QW = VY ^f f T^"*? d*. (5-22) 

ir exp L Jo 1 — o cos 

This formula can also be obtained as a particular case of the more general 
formula (2.24) by setting = 2ic in the upper limit of integration, utilizing 
(2.25), and changing the variable of integration by the substitution = 

0/2. 

Two partial checks orr any general formula for Q(R) = Q(R;b) or for 
Q*(R) = Q*{R;b) can be applied by setting b = and b =• 1, and comparing 
the resulting particular formulas with those obtained by integrating the 
formulas for P(R;0) and P'{R;l) obtained in Section 3, namely formulas 
(3.20) and (3.21) there. It is thus found that 

Q*(R;0) = exp(- J R 2 ) = [ P(R;0)dR, (5.23) 

J R 

Q(Ri ° = 2 fe f ^ [" I] "} = I" 1 "*' ° * (5 - 24) 
It will be recalled that the quantity between braces in (5.24) is extensively 
tabulated, and that it is sometimes called the 'normal probability integral.' 

Several of the above general formulas for Q(R) m p(R' <R) and for 
Q*{R) = p(R'>R) are closely connected with my 1933 paper. 27 Indeed, 
formulas (5.11), (5.14), (5.16), (5.19) and (5.22) above are the same as 
(53- A), (56- A), (52-A), (55-A) and (22-A), respectively, of the unpublished 
Appendix A to the 1933 paper; and (5.12), (5.13), (5.15), (5.17), (5.18) and 
(5.20) above were derived in the same connection, although they were not 
included in the Appendix A. 

Formula (5.22) was employed in the unpublished Appendix A of the 1933 
paper, being (22-A) there, as a basis for deriving two very different kinds 
of series type formulas for computing the values of p(R'>R) = Q*(R) 
underlying the values of p b , {R'>R) constituting Table I (facing Fig. 8) 
in that paper. 28 

26 This formula, (5.22), was derived by me in a somewhat different manner in the un- 
published Appendix A to my 1933 paper. Later I found that an equivalent formula, 
easily transformable into (5.22), had been given by Bravais as formula (51) in his classical 
paper "Analyse math6matique sur les probability des erreurs de situation d'un point," 
published in Memoires de l'Acad6mie Rovale des Sciences de l'lnstitut de France, 2nd 
series, vol. IX, 1846, pp. 255-332. (This is available in the Public Library of New York 
City, for instance.) 

« There the abbreviated symbols p(K' < R) and p(R' > R) were used with the same 
meanings as the complete symbols on the right sides of equations (5.1) and (5.2), respec- 
tively, of the present paper. 

28 Each of the two kinds of series type formulas comprised a finite portion of a con- 
vergent series plus an exact remainder term consisting of a definite integral. In the 
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In the present paper, formulas (5.11) and (5.16) have been used for numer- 
ical evaluation of Q(R) m p(R'<R) and of Q*(R) = p(R'>R) by numerical 
integration (employing 'Simpson's one-third rule'), aided by some of the 
considerations set forth in Appendix D. However, only a moderate number 
of values of these quantities have been thus evaluated — merely enough to 
afford a fairly comprehensive check on Table I of my 1933 paper, by means 
of a sample consisting of 60 values (about 26%) distributed in a somewhat 
representative manner over that table. These new values of Q*(R) = 
p(R'>R) = 1 - Q(R) are presented in Table 5.1 (at the end of this section) 
in such a way as to facilitate comparison with the old values, namely those 
in the 1933 paper. Thus, for any fixed value of R in Table 5.1, there are 
two horizontal rows of computed values of Q*(R), the first row (top row) 
coming from the 1933 paper, and the second row coming from the present 
paper. The third row of each set of four rows gives the deviations of the 
second row from the first row; and the fourth row expresses these deviations 
as percentages of the values in the first row. 

In the first row of any set of four rows, any value represents Q*(R) = 
p b (R'>R) obtained, in accordance with Eq. (22) of my 1933 paper, by 
adding exp (— R 2 ) to pb,o(R'>R) given in Table I there. In the second 
row of a set, any value represents Q*(R) = 1 — Q(R) as computed by for- 
mula (5.11) or (5.16) of the present paper: more specifically, the values for 
R = 0.2, 0.4, 0.6 and 0.8 were computed by (5.11); and the values for 
R =■ 1.6 and R = 2 by (5.16), taking r = 1/1.6 = 0.625 and r = 1/2 = 0.5 
respectively. 29 

In the 1933 paper, the values of p b {R'>R) = Q*(R;b) for b = and for 
b = 1 were omitted as being unnecessary there because their values could 
be easily obtained from the simple exact formulas to which the general 
formulas there reduced, for b = and 6=1. Those reduced formulas 
were the same as (5.23) and (5.24) of the present paper, except that (5.24) 
gives Q(R;l) instead of giving Q*{R;l) = 1 - Q{R;l). The values obtained 
from these two formulas, exact to the number of significant figures here 
retained, are given in Table 5.1 at the intersections of the first row of each 
set of four rows with the columns b = and b = 1. Therefore in these two 
columns the deviations (in the third row of each set of four rows) are devia- 
tions from exact values; the values in the second row of each set are, as 

use of such a formula for numerical computations, the expansion producing the con- 
vergent series was carried far enough to insure that the remainder definite integral would 
be relatively small, though usually not negligible; and then this remainder definite integral 
was evaluated sufficiently accurately by numerical integration. 

29 In the work of numerical integration, 'Simpson's one-third rule' was employed for 
R = 0.2, 0.4, 0.6, 0.8 and 2. For R = 1.6, so that r = 1/1.6 = 0.625, 'Simpson's one- 
third rule' was employed up to r = 0.620, and the ' trapezoidal rule' from r = 0.620 to 
r = 0.625. 
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already stated, those obtained by the methods of the present paper, employ- 
ing numerical integration. 

From detailed inspection of Table 5.1 it will presumably be considered 
that the agreement between the two sets of values of Q*(R;b) = p b {R'>R) 
is to be regarded as satisfactory, at least from the practical viewpoint, the 
largest deviation being less than one per cent (for R = 0.8, b = 0.9). 

Table 5.1 
Values of Q*(R) - p(R' > R) 



0.2 



0.4 



0.6 



0.8 



1.6 



2.0 






0.3 


.9608 


.9590 


.9623 


.9605 


.0015 


.0015 


.16 


.16 


.8521 


.8462 


.8537 


.8477 


.0016 


.0015 


.19 


.18 


.6977 


.6880 


.6992 


.6892 


.0015 


.0012 


.22 


.17 


.5273 


.5167 


.5290 


.5183 


.0017 


.0016 


.32 


.31 


.07730 


.07986 


.07727 


.07988 


-.00003 


.00002 


-.04 


.03 


.01832 


.02153 


.01823 


.02145 


-.00009 


- .00008 


-.49 


-.37 



0.4 



.9574 
.9590 
.0016 
.17 

.8410 
.8427 
.0017 
.20 

.6799 
.6814 
.0015 
.22 

.5081 
.5099 
.0018 
.35 

.08207 
.08210 
.00003 
.04 

.02394 
.02383 
.00011 
.46 



0.5 


0.6 


0.7 


0.8 


0.9 


0.95 


9550 


.9516 


.9463 


.9372 


.9168 


.8930 


9567 


.9528 


.9473 


.9387 


.9206 


.8925 


0017 


.0012 


.0010 


.0015 


.0038 


-.0005 


18 


.13 


.11 


.16 


.41 


-.06 


8335 


. 8228 


.8071 


.7830 


.7420 


.7127 


8351 


8240 


.8081 


.7841 


.7459 


.7125 


0016 


.0012 


.0010 


.0011 


.0039 


-.0002 


19 


.15 


.12 


.14 


.53 


-.03 


6686 


.6531 


.6324 


.6055 


.5721 


.5578 


6698 


.6540 


.6334 


.6065 


.5764 


.5572 


0012 


.0009 


.0010 


.0010 


.0043 


-.0006 


18 


.14 


.16 


.17 


.75 


-.11 


4969 


.4826 


.4656 


.4477 


.4316 


.4261 


4982 


.4840 


.4672 


.4488 


.4357 


.4266 


0013 


.0014 


.0016 


.0011 


.0041 


.0005 


26 


.29 


.34 


.25 


.95 


.12 


08522 


.0891 


.0938 


.0990 


.1042 


.1070 


08536 


.0892 


.0938 


.0989 


.1042 


.1069 


00014 


.0001 


.0000 


-.0001 


.0000 


-.0001 


16 


.11 


.00 


-.10 


.00 


-.09 


07681 


.0301 


.0337 


.0375 


.0414 


.0435 


02685 


.0302 


.0338 


.0376 


.0415 


.0436 


00004 


.0001 


.0001 


.0001 


.0001 


.0001 


15 


.3i 


.30 


.27 


.24 


.23 



1.00 



.84148 
.84124 
.00024 
.03 

.68916 
. 68897 
.00019 
.03 

.54851 
.54831 
.00020 
.04 

.42371 
.42355 
.00016 
.04 

.10960 
. 10958 
.00002 
.02 

.04550 
.04552 
.00002 
.04 



6. The Distribution Function For The Angle 

The distribution function P(0 | Ri 2 ) = P(0) for the angle of any complex 
variate W = R(cos 6+ ism 0) is denned by equation (1.10) on setting 
p = e, a = R, a = Ri = and <r 2 = R*. = °° J thus 

P(0)dd = p{6<6'<d+dd,Q<R'<™). (6.1) 

An integral formula for P(0) is immediately obtainable from (1.6) by 
setting p = 0, c = R, <r 3 = <n = R\ = and <r 4 = <r 2 = R» = °° 5 thus 



P(6) = [ P(R, 0) dR. 
Jo 



(6.2) 
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The rest of this section deals with the case where W = R(cos + i sin 0) 
is 'normal.' Since this case depends on b as a parameter, P(0) is here an 
abbreviation for P(d;b). 

A formula for P(6;b) m P(0) can be obtained by substituting P(R,6) 
from (2.15) into (6.2) and executing the indicated integration, which can 
be easily accomplished. The resulting formula is found to be 

^■sS- (6 - 3) 

This formula can also be obtained as a particular case of either of the 
more general formulas (2.19) and (2.20) by setting R = °o in (2.19) or 
R = in (2.20); also by adding (2.19) to (2.20) and then utilizing (1.10). 

In P(6) m P(d;b) it will evidently suffice to deal with values of in the 
first quadrant, because of symmetry of the scatter diagram. 

In P(6;b) it will suffice to deal with only positive values of b, as (6.3) 
shows that changing b to —b has the same effect as changing 20 to7r±20, 
or to ir/2±0; that is, P(6;-b) = P(tt/2±0;&). 

Fig. 6.1 gives curves of P(d;b), computed from (6.3), as function of 
with b as parameter, for the ranges 30 0^0^90° and O^b^l. 

The curves in Fig. 6.1 indicate that P(6;b) is a maximum at = 0° and 
a minimum at = 90°. These indications are verified by formula (6.3), 
as this formula shows that: 

Max P{d;b) = P(0°;b) = JL /^\^\ , (6-4) 

Min P(d;b) = P(90°;6) = 1 zj/pF? • (6-5) 



+ b 



Thence 



Min P(0;6)/Max P(6;b) = (l-6)/(l+6), (6.6) 

P(d;b)/MaxP{d;b) = P(d;b)/P(0°;b) = (l-b)/(l-b cos 20). (6.7) 

The curves in Fig. 6.1 indicate also that P(6;b) is independent of in 
the case b = 0. This is verified by formula (6.3), as this formula shows that 

P(B-fl) = 1/2tt. (6.8) 

Thence (6.3) can be written 

P(d- } b)/P(8;0) = ( V / r^T 2 )/(l-&cos20). (6.9) 

30 Beginning here, 6 will usually be expressed in degrees instead of radians, for prac- 
tical convenience. 
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By setting cos 20 = in (6.3), so that 6 = 45°, it is found that 

(\/r^T 2 )/2T = P(45°;b), (6. 



10) 



4? — <\jpWm<or~«o a> S 
odo'ddcidcid d d 




DISTRIBUTION FUNCTION, P(9;t>) 

Fig. 6.1— Distribution function for the angle. 

whence (6.3) can be written 

P(0;b)/P(4S°;b) = 1/(1-6 cos 20). 



(6.11) 
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REDUCED DISTRIBUTION FUNCTION, P(9; t>)/MAX P(9;D) 

Fig. 6.2 — Reduced distribution function for the angle. 
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In the case 6=1, the curves in Fig. 6.1 suggest, by limiting considera- 
tions, that P(0;1) is zero for all except = 0°, and that P(9;l) is infinite 
for = 0°. These conclusions are verified by formula (6.3), as this formula 
shows that: 

P(0;1) = for O°<0<18O°; P(0;1) = °o for = 0°, 180°. 

The curves in Fig. 6.1, though having the advantage of directly rep- 
resenting P(6;b) as function of with b as parameter, are somewhat trouble- 
some to use because of their numerous crossings of each other. This 
difficulty is not present in Fig. 6.2, which gives curves of P(0;6)/Max 
P(6;b), obtained by dividing the ordinates P(d;b) of the curves in Fig. 6.1 
by the respective maximum ordinates of those curves, as given by (6.4), 
so that the equation of the curves in Fig. 6.2 is formula (6.7). 

7. The Cumulative Distribution Function for the Angle 

The cumulative distribution function Q(<6,Ru) = Q{d) for the angle 
of any complex variate W = R(cos + i sin 0) is defined by equation 
(1.11) on setting p = 0, a = R, pi = 0i = 0, <n = Ri = and a 2 = Rz = °° ; 
thus 

Q(e) = p(o<e'<d, o<R'< oo). (7.1) 

A 'double integral' for Q{9), in the form of two 'repeated integrals,' can 
be written down directly by inspection of the p( ) expression in (7.1) 
or by specialization of (1.8); thus 



Q{e) = f 9 If P(R, 0) dRj dQ = f l^ P(R, 0) dd\ dR. (7.2) 
can be written formally as two 'single integrals,' 
Q(6) = [ P{6) dd = [ P(R\ < 0) dR, (7.3) 

f •'0 •'o 



Evidently these can be written formally as two 'single integrals,' 

I 

I Jo 

by means of the distribution functions P(0) = P(d\R n ) and P(R\ <0) 
given by the formulas 

P(0) = ( P(R, 0) dR, (7.4) P(R | < 0) = f P(R, 0) dd. (7.5) 

Jo J ° 

(7.4) is the same as (6.2). (7.5) is a special case of (1.6), and the left side 
of (7.5) is a special case of P(p \ <<r) defined by (1.13). 

The rest of this section deals with the case where W = R(cos + i sin 0) 
is 'normal.' Since this case depends on b as a parameter, Q{6) is here an 
abbreviation for Q(6;b). 

A natural and convenient way for deriving formulas for Q(d) is afforded 
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by the general formula (7.3) together with the auxiliary general formulas 
(7.4) and (7.5), beginning with the two latter. 

It will be convenient to dispose of (7.5) before dealing with (7.4), as (7.5) 
turns out to be the less useful. For when P(R,9) given by (2.16) is sub- 
stituted into (7.5), the indicated integration cannot be executed in general, 
as (7.5) becomes (2.18), wherin the indicated integration can be executed 
only for certain special values of the integration limit — by means of the 
special Bessel function formula (3.3). 

When P(R,6) given by (2.15), which is equivalent to (2.16) used above, 
is substituted into (7.4), it is found that the indicated integration can be 
executed, giving the previously obtained formula (6.3) for P(0) = P(0;b). 

A 0-integral formula for Q(9) = Q(6;b) can be obtained by substituting 
P{6) m P{d;b) from (6.3) into the first integral in (7.3), giving 

0(6- b) = Vlr=7T2 f dd = VJZI 2 f ** (76) 

VK ' 2tt Jo 1 - b cos 20 4tt J 1 - b cos * * " ' 

This formula can also be obtained as a particular case of the more general 
formulas (2.22) and (2.24) by setting R = °o in (2.22) or R = in (2.24); 
also by adding (2.22) to (2.24) and then utilizing (1.11). 

The integral in (7.6) is of well-known form, and the indicated integration 
can be executed, yielding the following two equivalent formulas for Q(8;b): 



Q(d; b) = i- 

27T 



tan" 



4tt 



— i 
cos 



(7.7) 
cos 26 — b 



Jl- 



.1 - b cos 20. 

In Q(6;b) it will evidently suffice to deal with values of in the first quad- 
rant, because of symmetry of the scatter diagram, and the resulting fact 
that Q(n 90°) = n/4, where ft = 1, 2, 3 or 4. 

In Q{9;b) it will suffice to deal with positive values of b, as (7.7) shows 
that 31 



Q(d; -b) = 



i-e(i±M)|. 



Fig. 7.1 gives curves of Q(6;b) = Q(6) computed from (7.7), as function 
of with b as parameter, for the ranges 0^0^90° and 0^6^ 1. 

Consideration of the scatter diagram of W or of its equiprobability curves, 
which are concentric similar ellipses, affords several partial checks on the 
curves in Fig. 7.1 and on formula (7.7) from which they were plotted. 

31 This relation can also be derived geometrically from the fact that the scatter dia- 
gram for —b is obtainable by merely rotating that for b through 90°, as shown by (2.6), 
or (2.7) and (2.8), or (2.11). 
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CUMULATIVE DISTRIBUTION FUNCTION, Q(9;b) 



Fig. 7.1— Cumulative distribution function for the angle. 

Thus, the fact that the curve for b = is a straight line, whose equation is 
Q(0;O) = 6/2tt = O /36O°, {b = 0), 
corresponds to the fact that for b = the equiprobability curves are circles. 
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The fact that the curve f or b = 1 is the straight line <3(0;1) = 1/4 = 0.25 
corresponds to the fact that for b = 1 the scatter diagram has degenerated 
to be merely a straight line coinciding with the real axis, so that no point 
outside of this line makes any contribution to Q(8;l). 

The fact that, at = 90°, Q(6;b) = Q(90°;b) has for all b the value 1/4 = 
0.25 corresponds to the fact that the area of a quadrant of the scatter 
diagram is one-fourth the area of the entire scatter diagram. Hence 
(3(360° ;b) - 4Q(90 o ;b) = 1, which is evidently correct. 
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APPENDIX A 

Derivation of Formula (2.15) for P{R,6) 

(2.15) will here be derived from (2.11) by utilizing the fact that the 'areal 
probability density', G, at any fixed point in the scatter diagram must be 
independent of the system of coordinates; for G dA gives the probability 
of falling in any differential element of area dA, and this probability must 
evidently be independent of the shape of dA (assuming that all linear dimen- 
sions of dA are differential, of course). Thus, indicating the element of 
area by an underline, we have, in rectangular coordinates, 

G dUdV = P(U,V)dUdV, (Al) whence G = P(U,V). (A2) 

In polar coordinates, 

G RdddR = P(R,e)dRd6, (A3) whence G = P(R,d)/R. (A4) 

Comparing these two expressions for G shows that 

P(R,6) = RP(U,V). (A5) 

Thus, a formula for P(R,6) can be obtained from (2.11) by merelv multiply- 
ing both sides of that formula by R. However, in the resulting formula it 
will remain to express V and T' in terms of R and 0, by means of the relations 

U = R cos 0, (A6) V = R sin 0. (A7) 

The final result, after a simple reduction, is (2.15), which is thus proved. 

APPENDIX B 

Formulas of the Curves in Fig. 3.3 

As in equation (3.22), R c will here denote the critical value of R, that is, 

the value of R at which P(R) = P(R;b) has its maximum value; and T c 

32 Formula (A5) can be easily verified by the entirely different method which utilizes 
(1.23). 
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will denote the corresponding value of T, whence T c is given in terms of 
R c and b by (3.22). 

A formula for dP{R)/dR could of course be obtained directly from (3.4) 
but it will be found preferable to obtain it indirectly from the less cumber- 
some formula (3.8) containing the auxiliary variable T denned by (3.6). 
Evidently, since b does not depend on R, 



dP{R) = dP(R) dT = 2bR dP(R) 
~~dR dT dR 1 - b 2 dT 



(Bl) 



Thus, since the factor 2bR/(l-b 2 ) cannot vanish for any value of R (except 
R = 0), the only critical value of R must be that corresponding to the value 
of T at which dP(R)/dT vanishes, namely T c , since T c has been defined 
to be the value of T corresponding to R c . (Incidentally, equation (Bl) 
shows that T c is equal to the value of T at which P(R) is an extremum 
when P(R) is regarded as a function of T.) From (3.22), 



0*2) 



Rc T c 



1- b* b 

Evidently T c and R c must ultimately be functions of only b. The next 
paragraph deals with T c , which evidently has to be known before R c can 
be evaluated. 

From (3.8) it is found that, since dI {T)/dT = Ii(T), 

*m = P( R )\± + !*n _ n . (B3) 

dT K } 12T T h(T) b\ 

Hence, since P(R) does not vanish for any value of R (except R = and 
R = oo ), T c will be a root of the conditional equation obtained by equating 
to zero the expression in brackets in (B3). This conditional equation is 
transcendental in T c and apparently has no closed form of explicit solution 
for T c ; and its solution by successive approximation, or otherwise, would 
likely be rather slow and laborious. However, the bracket expression in 
(B3) shows that b can be immediately expressed explicitly in terms of T c 
by the equation 

b = ?I2 . (B4) 

1 + 2T c h(Tc)/h{Tc) 

For some purposes, the following two equations, each equivalent to (B4), 
will be found more convenient: 

T = 2 + Tc u^y (E5) 

l£ = Ml . (B6) 

b 1 - bh(T c )/h{T c ) 
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On account of (B2), the right sides of (B5) and (B6) are equal not only to 
T c /b but also to i? 2 c /(l-* 2 )- 

Since the utilization of formulas (B4), (B5) and (B6) for computing the 
curves in Fig. 3.3 will involve taking T c as the independent variable and 
assigning to it a set of chosen numerical values, the natural first step is to 
find approximately the range of T c corresponding to the 6-range, OSa&^l, 
in order to be able to choose only useful values of T c . This step will be 
taken in the next paragraph. 

Equation (B6) shows that T c /b = 1/2 when b = 0, and hence that T c = 
when b = 0; and this last is verified by (B4). The other end-value of the 
r c -range, namely the value of T c for b = 1, cannot be found explicitly 
and exactly. However, rough values of limits between which it must lie 
can be found fairly easily as follows: To begin with, each of the equations 
(B5) and (B6) shows that T c t b/2, for all values of b in 0=6^ 1; in par- 
ticular, T c > 1/2 when 6=1. An upper limit for T c for any value of 
b can be found from (B5) by utilizing the power series expressions for 
I\{Tc) and Io(T c ), whereby it is found that 

£S4 =JI~,' (B7) where H = 1 - % c < 1. (B8) 

io{l c) £ o 

On substituting (B7) into (B5) and then solving for T c in terms of b and 
H, it is found that 

To = 6/(1 + Vl - Hb>). (B9) 

On account of (B8), (B9) shows that 

T c < 6/(1 + Vr^i" 2 ), (BIO) 

whence, in particular, T c <\ when 6=1. By successive approximation 
or otherwise, it can now be rather quickly found that, when 6 = 1, To — 
0.79 (to two significant figures). 33 

From the preceding paragraph, it is seen that, when 6 ranges from to 1, 
T c ranges from to about 0.79; T c /b ranges from 0.5 to about 0.79; and, 
on account of (B2), R c ranges from \/0.5 = 0.707 down to 0. 

The curves in Fig. 3.3 are constructed with the aid of the formulas and 
methods of this appendix as follows: First, a set of values of T c is chosen, 
ranging from to 0.79 and slightly larger. Second, for each such chosen 
T c the right side of (B5) is computed, thereby evaluating T c /b and also 
R c /{l — b 2 ), these two quantities being equal by (B2). Third, the cor- 
responding value of 6 is found by dividing T c by T c /b; less easily, it could 

33 Because of the special importance of b = 1 in other connections, T c for b = 1 was 
later evaluated to four significant figures and found to be T c = 0.7900; thence, by sub- 
stituting this value of T into (3.8), along with b = 1, it was found that Max. P(R;l) 
= 0.9376, which occurs at R = R c = 0, by (B2). 
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be found by substituting T c into (B4). Fourth, fr om T c /b the value of 
\/T c /b is found, and thereby the value of R c /Vl — b 2 and thence R c . 
Finally, Max. P(R;b) is computed by inserting the critical values into any 
of the various (equivalent) formulas for P(R;b), namely (3.4), (3.7), (3.8), 
(3.10) or (3.12). 

APPENDIX C 

FOMULAS OF THE CURVES IN FlG. 4.3 

The first six equations of this appendix are given without derivation 
and almost without any comments because they correspond exactly and 
simply to the first six equations, respectively, of Appendix B. Beginning 
with the second paragraph of the present appendix, the close correspondence 
ceases. 

dP(r) _ dP{r) dT = -26 dP(r) 

~dr dT dr " = (1 - 6*)r» dT 

1 T e 



(CI) 
(C2) 



(1 - b 2 )rl ~ b ■ 
dP(r) .. pM f 3 , h(T) _ 1] (C3) 

h _ ^Tc /p^\ 

3 + 2T C h{T c )/h(T e ) ' 

I± = l + Tc h{ - Te) . (C5) 

b 2 + ° h(T e ) 

I± = 3 / 2 . (C6) 

b 1 - bIi(T e )/I (T e ) 

The bracketed expression in (C3) is seen to be obtainable from that in (B3) 
by merely changing T to T/3 wherever T does not occur as the argument 
of a function; hence the three equations following (C3) are obtainable from 
the three equations following (B3) by correspondingly changing T c to 
T c /3. (In this appendix, as in Section 4, small c is purposely used as a . 
subscript to indicate a 'critical' value, whereas in Section 3 and in Appendix 
B, capital C is used for that purpose.) 
For use below, it will here be noted that 

h{T c )/h(T c ) = Ni{T e )/No{T e ), (C7) 

as will be seen by dividing (3.16) by (3.13). On account of (3.17) and (3.14), 
(C7) shows that for large values of T c the right side of (C7) is equal to 1 
as a first approximation, and to 1 - 1/2 T c as a second approximation; 
thus, for large T e , 
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h{T e )/h(T c ) * 1 - 1/2T C = 1. (C8) 

The first step toward computing the curves in Fig. 4.3 is to find approxi- 
mately the TV range corresponding to the 6-range, O^b^l. This is done 
in the course of the next four paragraphs. 

When b = 0, equation (C6) shows that T c /b = 3/2 and hence thai 
T c = 0; or, what is equivalent, b/T c = 2/3 and hence 1/T C = co (since 
b= 0). 

When b = 1, T e = oo, as can be easily verified from equation (C4), 
(C5) or (C6) by utilizing (C8). 

Thus, from the two preceding paragraphs, it is seen that, when b ranges 
from to 1, b/T c ranges from 2/3 to 0; T c /b from 3/2 to <x> ; and T c from 
to co. 

Since T c = <x when b = 1, the choosing of a set of finite values of T c 
will necessitate an approximate formula for computing T e for values of 
b nearly equal to 1, which means for very large values of T. Such a formula 
is easily obtainable from (C5) by utilizing the approximation 1 — 1/2 T c 
in (C8), whereby it is found that, for large T c , 

T e = b/(l-b), (C9) b/T e » l-b. (CIO) 

As examples, these approximate formulas give: When b = 0.99, T c = 99, 
b/T c = 0.01; when b = 0.9, T c = 9, b/T c = 0.1. It will be found that 
even in the second example the results are pretty good approximations. 

The curves in Fig. 4.3 are constructed with the aid of the formulas and 
methods of this appendix as follows: First, a set of values of T e is chosen, 
ranging from to about 100 (the latter figure corresponding approximately 
to b = 0.99). Second, for each such chosen T c the right side of (C5) is 
computed, thereby evaluating T c /b and also 1/(1 — b-)r 2 c , these two quan- 
tities being equal by (C2). Third, the corresponding value of b is found 
by dividing T c by T c /b; less easily, it could be found by substituting T c 
into (C4). Fourth, from T c /b the value of y/T c /b is found, and thereby 
the value of \/r c s/l — b 2 and thence r c . Finally, Max P(r;b) is computed 
by inserting the critical values into any of the (equivalent) formulas for 
P(r;b), namely (4.2), (4.3) or (4.4). 

APPENDIX D 

Some Simple General Considerations Regarding the Evaluation of 
Cumulative Distribution Functions by Numerical Integration 

This appendix gives some simple general considerations and relations 
that may sometimes facilitate and render more accurate the evaluation 
of cumulative distribution functions by numerical integration. 
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Some of these considerations and relations have found application in 
Section 5 in the evaluation of the cumulative distribution function for the 
modulus R = | W |. For this reason, the variate in the present section 
will be denoted by R, though without thereby restricting R to denote the 
modulus; rather, R will here denote any positive real variate, though it 
should preferably be a 'reduced' variate, so as to be dimensionless, as in 
equation (2.9). The restriction of R to positive values is imposed because 
it is strongly conducive to simplicity and brevity of treatment, without 
constituting an ultimate limitation. The reciprocal of jR will be denoted 
by r, as previously. 34 

We may wish to evaluate numerically the cumulative distribution func- 
tion p(R'<R) = Q(R) or p(R'>R) = Q*{R) or both. Since these are not 
independent, their sum being equal to unity, the evaluation of either one 
determines the other, theoretically. However, when the evaluated one is 
nearly equal to unity, the remaining one may perhaps not be evaluable 
with sufficient accuracy (percentagewise) by subtracting the evaluated one 
from unity. Then it would presumably be advantageous to introduce 
for auxiliary purposes the variable r = l/R, since evidently 

p(R'>R) = p(l/R'<l/R) = pir'Kr), (Dl) 

p(R'<R) = p(r'>r) - 1 - p(/<r). (D2) 

Thus, if p{R'>R), in (Dl), is small compared to unity, it is presumably 
evaluable with higher accuracy percentagewise by dealing with p(r' <r) 
than with 1 - p{R'<R). Incidentally, after p(r'<r) has been evaluated, 
it might be used in (D2) to arrive at a still more accurate value of p(R'<R) 
than had originally been obtained directly by numerical integration. 

Assuming that we have a plot (or a table) of the distribution function 
P(R), we can evidently evaluate 

P{R'<R°) = / P(R)dR (D3) 

Jq 

directly by numerical integration, provided the plot is sufficiently extensive 
to include R°; if not, we can, by (D2), resort to 



p(R'<R°) = 1 - p(r'<r°) = 1 - / P(r)dr, 

Jo 



(D4) 



assuming that a sufficiently extensive plot (or table) of P(r) is available 
and applying numerical integration to it. 
Even if the plot of P(R) used in (D3) is sufficiently extensive to include 

34 The restriction of R, and hence of r, to positive values is seen to be absent from equa- 
tions (Dl), (D2), (D5) and (D6) but present in (D3), (D4), (D7) and (D8). 
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R , so that (D3) could be evaluated, it might be that (D4) would result 
in greater accuracy; this would presumably be the case when p{R'<R ) 
is nearly equal to unity. 
Evidently an evaluation of 

p(R'>R°) = [ P(R)dR (D5) 

directly by numerical integration would be less satisfactory than the evalua- 
tion of p(R'<R°) in the preceding paragraph. For, due to the presence 
of the infinite limit in the integral in (D5), the plot of P(R) would have to 
be carried to a large enough value of R so that the integral from there to oo 
would be known to be negligible. This difficulty can be avoided by start- 
ing with the relation 

p(R'>R°) = 1 - p{R'<R°) (D6) 

and substituting therein the value of p{R'<R ) given by (D3) or (D4), 
resulting respectively in the following two formulas: 



r R0 
p(R'>R°) = 1 - P(R)dR, 
Jo 



(D7) 



p(R'>R°) = p(r'<r°) = f P(r)dr, (D8) 

jo 

the integrals in which are evidently suitable for evaluation by numerical 
integration, none of the integration limits being infinite. If p(R'>R°) 
is small compared to unity, (D8) would presumably be more accurate 
(percentagewise) than (D7). If the plot of P(R) is not sufficiently exten- 
sive to include R°, (D7) evidently could not be used; but, instead, (D8) 
could be used if the plot of P(r) were sufficiently extensive to include r . 
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